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Abstract 

Experimental results for the reflection coefficient of shock-compressed dense Xenon plas- 
mas at pressures of 1.6 - 17 GPa and temperatures around 30 000 K using a laser beam 
with A = 1.06 fim are compared with calculations based on different theoretical approaches 
to the dynamical collision frequency. It is found that a reasonable description can be given 
assuming a spatial electron density profile corresponding to a finite width of the shock wave 
front of about 2 • 10^^ m. 

1 Reflectivity measurements 

Experiments with an explosively driven generator of shock waves to produce dense 
nonideal Xenon plasmas were reported in |T|. The reflectivity was measured using a 
laser system with the wave length = 1.06 /im. The results of the experiments are 
shown in Tab. 1. The thermodynamic parameters of the plasma were determined from 
the measured shock wave velocity. The plasma composition was calculated within a 
chemical picture 0. Working with a grand canonical ensemble virial corrections 
have been taken into account due to charge-charge interactions (Debye approxima- 
tion). Short range repulsion of atoms and ions (including multiply charged positive 
ions) was considered via second and third virial coefficients within a virial expansions. 
In the parameter range of the shock wave experiments, derivations of up to 20 % for 
the composition have been obtained depending on the approximations for the equation 
of state. This is within the accuracy of the experimental values of the reflectivity. 
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0.26 


10.5 


29250 


2.70 


7.1x10^1 


5.4x10^1 


0.36 


16.7 


28810 


3.84 


9.1x10^1 


8.6x10^1 


0.47 



Tab. 1: Experimental results for the reflectivity R of Xenon plasmas at given parameter 
values: pressure P, temperature T, mass density p, free electron number density Ue and 
density of neutral atoms n^. 
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It was suggested that the measurement of reflectivity allows to determine the con- 
centration of free electrons in the plasma. However, investigating these data 0] , it was 
not possible to find a direct relationship between the values of free electron density 
He and the reflectivity R. The reflection coefficient increases smoothly with the free 
electron density. It approaches only slowly values characteristic for metals, although 
the critical density for metallic behaviour rf^ = 1.02 x lO^^cm"'^, where the plasma fre- 
quency, = ^Uec'^/ {eonie) with me the electron mass, coincides with the frequency 
uji of the probing laser pulse, is exceeded even at the lowest density. However, this 
critical density is taken from the RPA approximation for the dielectric function in the 
long-wavelength limit e = 1 — u^i/ uf where total reflection occurs due to a vanishing 
dielectric function. Taking collisions into account this will be smoothed out. On the 
other hand, a sharp boundary between plasma and undisturbed gas in front of shock 
wave is assumed. 

In Jl]], the spatial structure of the ionizing shock wave was discussed, showing three 
characteristic zones, which may influence the electromagnetic wave propagation. In a 
precursor zone, the gas is heated, but the influence on the wave propagation is small. 
In the region of the shock wave front a steep increase of the free electron density is 
expected. The width of the wave front is determined by relaxation processes in the 
plasma. It was estimated to be of the order d ~ 10"'' m [|l|, which is one order of 
magnitude less than the laser wave length. Under this condition rf/A -C 1, the laser 
beam reflection was assumed to be determined by the electron properties of the plasma 
behind the shock wave front. An expression derived from the Fresnel formula in the 
long-wavelength limit 



6(^) - 1 ' 



eM + 1 



(1) 



was applied where the frequency u has to be taken at the laser frequency ui = 1.8 ■ 10^^ 
Hz. The complex frequency-dependent dielectric permittivity 

e(c^) = 1 + —^{u;) = 1 - (2) 
eouJ uj[uj + ti'[uj)\ 

has been related to the dynamical conductivity cr{uj) or the dynamical collision fre- 
quency ^{uj). The Drude formula 

a(uj) = —— — - — (3 
^ ' z/(0) -iu ^ ^ 

follows if the collision frequency is taken in the static limit i/(0) = eouj'^i/a{0) relating 
it to the static conductivity o"(0). Different expressions for a{0) have been considered 
in |l[], but no satisfying explanation of the experimental results has been given there. 
As will be seen in the following discussion, the Drude model is not an appropriate 
approximation under the experimental conditions considered. In the present work, 
improvements in the calculation of the dielectric function as well as considerations 
concerning the shape of the wave front will be discussed to flnd a consistent theoretical 
approach to the measured reflectivities. 
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2 Reflection by a step-like plasma front 

In this section, we investigate the reflection coeflicient at a shock wave front where its 
width d shall be neglected. According to the Fresnel formula, the reflection coeflicient 
at a step-like plasma front @] 



R{uj) 



1 - Z{uj) 



Z{uj) 



(4) 



is related to the surface impedance Z (normal incidence). The surface impedance can 
be related to the transverse dielectric function via 

iuj 1 
^("^ = --C l-J^ k^-u-e,{k,u)/c- ■ 

This is a more general expression than (|l|) for the reflectivity. Assuming that the 
transverse dielectric function et{k,u) is independent of the wave vector k, the integral 
@ can be evaluated leading to the expression (|l]) given above. 

In general, the dielectric function is related to a dynamical, nonlocal collision fre- 
quency ^{k, u), 

et{k,uj) = 1 — J- , (6) 

uj[uj + zutyk, u!)\ 

which is deflned by extending the above deflnition to flnite values of the wave 
vector k. In recent papers 0,0, an approach to the dielectric function within a linear 
response theory was developed. Different approximations have been investigated with 
respect to their consistency. Having this in mind, we are able to give a more precise 
description of the dielectric function at the plasma parameters considered here. 

2.1 Born approximation 

Before considering the nonlocal dielectric function below, we discuss the long-wavelength 
limit k ^ 0, replacing et{k,u!) by et{0,uj) in Eq. @. In this case, the transverse and 
longitudinal dielectric function are identical. The dynamical collision frequency can 
be evaluated in Born approximation with respect to the statically screened potential 
(Debye potential), see as (the non-degenerate case is considered) 



Born / \ 
U [Ujj = —ig Tie 



dy _ / dxe~^''~y^ = — , (7) 

{n + y'^Y J -oo xy{xy — uj — irj) 



where 



n - 



8eome(/cijT)2 

9 



24v^7r5/2e2 



m, 



1/2 



e 



and u = Tioj / {AkBT). The second integral in (^ is a complex quantity, 

dxe-^^-y^"v^-^ — - + tne-^^/y-y)'^ — . 

D xy[xy — uj) yuj 
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Evaluating the collision frequency in the static limit {u = 0), the result is the 
Ziman formula applied to a Debye potential in the case of non-degeneracy, 

i^Born^Q) = 47r(7ne H dy ^ c-y' . (8) 

In the following Tab. 2, the resulting values B^°^'^ for the reflectivity calculated from 
the Drude formula (|^) with the static collision frequency (§) are compared with the 
reflectivities i^^^'^ resulting from the dynamical collision frequency in Born approxi- 
mation (0) taken at the laser frequency uji. 

Compared with the observed reflectivities, the results obtained in Born approxima- 
tion are too small. Similar values were also reported in As well known, the Born 
approximation (Faber-Ziman result) underestimates the value of the dc conductivity. 
The correct low-density value of a^c is given by the Spitzer formula and can be ob- 
tained using a renormalization factor as discussed in see also discussion in Sec. 2.3 
below. The use of the dynamical conductivity increases the reflectivity by about 15 
%, but it also fails to produce the steep dependence of the reflectivity on the electron 
density as observed in the experiment, see Tab.|l|. 



2.2 Nonlocal conductivity 

The general expression (^), (|^) for the reflection coefficient contains an integral over the 
k dependent transverse dielectric function. We will discuss the effect of the nonlocal 
conductivity using the Mermin approximation as given in ^. The RPA solution is 
extended by introducing a complex frequency argument which contains the collision 
frequency. The Mermin expression for the dielectric function obeys particle number 
conservation. For the transverse dielectric function we obtain [Bl 



^Mermin / 1 
-t l'^) 



Mermin 



{k,uj) 




1 - 



1 



Mermin 



{k,uj) 



(9) 



with 



^Mermin 



{k,uj) 



(l+z^) [er^{k,u; + tu{u;)) 



1 + [e^Hk, uj + zuioo)) - 1] / [e^^k, 0) - 1] ' 



P,GPa 


UJpl/uJl 


pBorn 
-"-de 


^Born 


pMermm,Born 


1.6 


1.33 


0.272 


0.304 


0.272 


3.1 


1.77 


0.342 


0.351 


0.351 


5.1 


2.10 


0.381 


0.380 


0.376 


7.3 


2.36 


0.404 


0.399 


0.409 


10.5 


2.64 


0.429 


0.419 


0.443 


16.7 


2.99 


0.457 


0.447 


0.478 



Tab. 2: Reflectivities from step-like density profiles, calculated at the experimental parameter 
values. R^c'^'^ - Drude formula (^) with static collision frequency (^), R^°^^ - dynamical 

collision frequency in Born approximation (^), R^^^"^^^' - Eq. (^) with the Mermin 
nonlocal dielectric function and static collision frequency in Born approximation. 



H. Reinholz, Plasma Reflectivity 



5 



1 + [fi^^'Hk, ^ + ii^iuj)) - 1] / [fi^P^ik, 0) - 1] ' 



/^^"^(fc,^) = r^i 



The RPA expressions for the dielectric function of a Maxwellian plasma are 

efP^(A;,cu) = 1 + ^ (2 + ZeD{ze) + z,D{z,)) , (10) 



A;2 



with 



and the Dawson integral 

D{z) = ivr^/^g-^^ [1 Erf (iz)] . (12) 

Using the static collision frequency (||), the nonlocal dielectric function was calcu- 
lated. The results ^j-q q\^q given in Tab. 2. There is no essential modifi- 
cation if the fc-dependence of the dielectric function is taken into account. This can 
be explained by the fact that, at the conditions given, the main contributions to the 
integral over k come from the region of ~ O.OOOla^j^. A comparison of the mean free 
path and the much larger skin depth also shows, that nonlocal effects are not relevant 
here. We conclude that neither the account of the dynamical properties of the collision 
frequency nor the account of the k dependence in the dielectric function will essen- 
tially modify the calculated reflectivity. In particular, the discrepancy between the 
calculated and measured reflectivity as a function of the electron density, see Tab.^ 
cannot be cured. 



2.3 Strong collisions and renormalization 

The Born approximation can and should be improved considering strong collisions and 
renormalization as described, e.g., in 0. Restricting ourselves to the static case, we 
discuss different approximations for a^c- We introduce dimensionless parameters F, 
characterizing the non-ideality and degeneracy, respectively, of the plasma: 



r 



Ai^t^ksT \ 3 



1/3 



e 



2mf.kBT 



(Svr^rig 



\-2/3 



The dc conductivity for the statically screened Coulomb potential in Born approx- 
imation (7^°'^'^ was already given in Sec. 2.1 via (^. In the low-density limit, the 
dc-conductivity reads 



^Born 



0.299 



{AneonksTf/' 



1/2 



2 r 



(13) 
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However, the correct expression for the dc conductivity in this hmit is given by the 
Spitzer formula 



Spitzer 
'"dc 



0.591 



g2^1/2 



- inr 

2 



(14) 



It is obtained by taking into account strong colhsions which modify the Coulomb 
logarithm in (|1^). The prefactor also changes from taking into account higher moments 
of the distribution function. 

Recently, an interpolation formula for the dc conductivity of a fully ionized Coulomb 
plasma was derived [H, 



cr, 



ERR 
dc 



03/2 



Dln{l + A + B) - C - 



62 + re 



(15) 



where T in K, a in (fim) ^, and with the functions 

3 1 + «4/r^e 
1 + a2/r2e + 03/1492 
63(1 + C3e)/re/(i + 036^/^), 
c4/(in(i + r-i) + c5r2e). 



A 

B 
C 
D 



[ai + ci ln(c2r=^/2 ^ ^) 



The set of parameters is given by ao = 0.03064, oi = 1.1590, 02 = 0.698, 03 = 
0.4876, 04 = 0.1748, 05 = 0.1, = 0.258, hi = 1.95, 62 = 2.88, 63 = 3.6, Ci = 1.5, 
C2 = 6.2, C3 = 0.3, C4 = 0.35, C5 = 0.1. They are fixed by the low-density expansion of 
the dc conductivity ([l^) , the strong degenerate limit and numerical data in for the dc 
conductivity the intermediate parameter region. 

Using the Drude formula (^ with the conductivities crf^^ we obtain the reflectiv- 
ities -Rd(?^5 see Tab. |. For comparison, experimental values for the dc conductivity 
[|T0[] are given. The fit formula seems to overestimate the conductivity. The results for 
the reflectivities R^^ are rather high, even exceeding the measured values. Obviously, 
the size of the calculated reflectivities can be shifted considerably according to the 
approximations made for the dc conductivity. But the strong increase of the measured 
values can not be explained yet. 



P/GPa 




Born 
^dc 


'^dc 


dEHH 

-"-dc 


1.6 


72 000 


45 300. 


90 000 


0.502 


3.1 


82 000 


59 400. 


125 000 


0.586 


5.1 


97 000 


72 400. 


160 000 


0.629 


7.3 


97 000 


83 100. 


195 000 


0.660 


10.5 


97 000 


95 500. 


240 000 


0.691 


16.7 


100 000 


114 000. 


311 000 


0.728 



Tab. 3: Values of the dc conductivity in (fim)"^ for different approximations. (Texp - 
experimental estimates [^], - static Born approximation (P), c^^?-^ - interpolation 

formula (|^) and corresponding reflectivity -Rf;?"^ calculated using Drude formula (|3|). 
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It is possible to extend the approach given here for the dc conductivity to the 
dynamic conductivity, as shown in ^ using the Gould-DeWitt approach. However, 
as already shown for the Born approximation, which is part of the Gould-DeWitt 
approach, we do not expect a significant modification of the density dependence of 
the reflection coefficient on the electron density. We also do not expect a significant 
modification if the nonlocal, k dependent dielectric function is considered, because 
only the long-wavelength limit contributes in calculating the impedance. Furthermore, 
within a more accurate treatment the contribution of the interaction with neutral 
atoms should be included as well. However, at the temperatures considered here the 
influence of neutrals on the conductivity is small. 

In conclusion of this Section, we could show that improvements in the theory of 
the dielectric function lead to substantial modifications in the refiectivity. For the 
parameter values considered here, the theory predicts a high value of the refiectivity 
close to the values obtained from the interpolation formula ([T5|) . However, we are still 
not able to describe the strong variation of the measured refiectivity with the electron 
density. 



Obviously, it is not possible to interpret the measured values of the refiection coeffi- 
cient of dense Xenon plasmas within the assumption of a step-like density profile, i.e. 
d ~ for width of the shock wave front. In particular, the steep increase of the re- 
fiectivity only at densities above the critical one can not be explained despite a highly 
sophisticated approach to the calculation of the dielectric function. 

Using the interpolation formula for the dc-conductivity ([T5|) and the Drude formula 
@, effective densities necessary to reproduce the measured values for the reflection 
coefficient can be deduced. They have been found to lie between 0.75 and 1.6 of the 
critical density = 1.02 x 10^^ cm~'^ where the plasma frequency coincides with the 
frequency of the probe laser. 

These effective densities can be considered as an argument that the reflection of 
electromagnetic radiation occurs already in the outer region where the density is low. 
Within a simplifled picture, considering a proflle with increasing electron density, the 
radiation penetrates the low-density region of the plasma up to the region where the 
density approaches the critical value. Here the wave will be reflected. 

To perform an exploratory calculation, a density proflle was assumed where the 
electron density increases linearly with distance z from zero up to the saturation value 
rie at the distance d. We assume the following linear dependence of the dielectric 
function on the distance of the shock wave front z 



L is the depth where the critical density is reached and the radiation is assumed 
to be reflected. The total width of the shock wave front is than determined by 



The collision frequency i^cr = 4.09 x 10^^ s ^ was determined from the static colli- 
sion frequency at the critical density using the interpolation formula for the static 



3 Smooth density profile 




(16) 



d 




(17) 
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P,GPa 


R 


L/(10-'m) 


ci/(10-''m) 


4n/(10-^m) 


1.6 


0.096 


6.01 


1.13 


1.26 


3.1 


0.12 


5.44 


1.82 


2.08 


5.1 


0.18 


4.40 


2.07 


2.54 


7.3 


0.26 


3.46 


2.06 


2.68 


10.5 


0.36 


2.62 


1.95 


2.56 


16.7 


0.47 


1.94 


1.84 


2.12 



Tab. 4: Calculation of the width d of the shock wave front from the reflectivity R, assuming 
a linear density profile. The critical density ricr is reached at depth L^r- dcr according to 
Eq. (|l^), diin solving the propagation of radiation in a linear density profile. 



conductivity (|T5|). The reflectivity can be calculated via ( see [pT[] 

8z/f.r-^' 



Taking the experimental values R for the reflectivity, the width dcr of the shock wave 
front according to (0) is given in Tab.H. 

Compared with the value given in , see also Sec. 1 above, our estimation of the 
width of the shock wave front is larger by about one order of magnitude and almost 
independent of the thermodynamic parameters. Only the value at the lowest density 
value comes out to be smaller. 

A more rigorous treatment of the propagation of laser radiation through a shock 
wave front should take into account the dependence of the collision frequency on the 
local electron density. An arbitrary density profile n^iz) can be approximated by a 
sufficiently large number of thin layers with constant electron density, and solving the 
boundary conditions when going from one slab to the next one. In such a way, any 
arbitrary density proffie ne{z) can be calculated. 

We have considered a linear dependence of the electron density on the distance z 
within the shock wave front of width d, ne{z) = Ugz/d, with being the electron den- 
sity behind the shock front. The subdivision of the width d into equidistant slabs has 
been increased until convergence was reached. For given density rig, the interpolation 
formula for the dc conductivity crf^^ has been used to find the dielectric function of 
the corresponding slab according to the Drude formula (^. Results d^^ reproducing 
the experimental reflectivity values R are shown in Tab.^ Convergence was reached 
by dividing dn^ into 16 equidistant slabs. 



4 Conclusion 

In order to infer plasma parameters from optical reflection coefficient measurements 
of dense Xenon plasmas, we propose a width of the shock wave front of about 2 ■ 10~^ 
m. Only at the lowest density measured, the value for d may be smaller. The linear 
density proflle should be considered as a model only to estimate the width of the 
shock wave front. Our method to approximate a given density proflle by a sufficiently 
large number of thin layers with constant electron density can be applied to any ne{z) 
which, in general, could be obtained by determining the distribution function for the 
non-equilibrium process of the shock wave propagation. 
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Our calculations are based on an interpolation formula for the dc conductivity, 
obtained from a systematic quantum statistical treatment of limiting cases. In partic- 
ular, the account of the renormalization factor and of strong collisions is essential to 
obtain the correct low-density limit. The uncertainty in using the interpolation for- 
mula increases for F < 1. The values for the plasma parameter F for the dense Xenon 
plasmas considered here are in the region 1 < F < 2, and 5 > O > 1.5 and we esti- 
mate the error of about 30 %. Furthermore, the dynamical collision frequency should 
be used instead of the static one. Using the Gould-DeWitt approach, the dynamical 
collision frequency was investigated in and minor modifications (about one half of 
that given above) are expected. The effects of non-locality can be neglected in the 
region of plasma parameters considered, as shown above for the Born approximation. 

Xenon under the conditions considered is a partially ionized plasma. The com- 
position is shown in Tab.|I]. The conductivity as well as the related quantities are 
influenced by the neutral component, which leads to a modification of the reflection 
coefficient. However, at the temperatures of about 30 000 K the contribution of the 
neutral component to the conductivity is small and will not be considered here. Simi- 
larly, the role of non-equilibrium effects such as relaxation of the composition will also 
not be considered here. 

Shock wave front investigations may be the subject of forthcoming experimental 
work. An interesting point would be the simultaneous determination of reflectivities 
at different frequencies. It is expected that in this way more information about the 
density profile of the shock wave front can be obtained. 
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